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Abstract. Given two point sets S and T, in a many-to-many matching between S and T each 
point in S is assigned to one or more points in T and vice versa. A generalization of the many-to- 
many matching problem is the limited capacity many-to-many matching problem, where the number 
of points that can be matched to each point, that is the capacity of each point, is limited. In this 
paper, we provide an O (« 2 ) time algorithm for the one dimensional minimum-cost limited capacity 
many-to-many matching problem, where |S| + \T\ = n. Our algorithm improves the best previous 
time complexity of 0(fcn 2 ), that in which k is the largest capacity of the points in S U T. In this 
problem, both S and T lie on the real line and the cost of matching s £ S to t 6 T is equal to the 
distance between s and t. 

Key words, many-to-many point matching, one dimensional point-matching, limited capacity 
point matching 

1. Introduction. Given two point sets S and T with |5| + \T\ = n, a matching 
function between S and T assigns some points of one set to some points of the other set 
PQ . The matching has applications in a variety of different fields such as computational 
biology [2], operations research [3], pattern recognition [4], computer vision [5]. The 
matching problem has many variations. In a one-to-one matching between two point 
sets, each point of one set is matched to a unique point of the other set, so there 
exists a perfect correspondence between the members of the sets [6]. In a many-to- 
one matching from S to T, each point of S is assigned to one or more points of T, and 
each point of T is assigned to a unique point of S [7\. In a many-to-many matching 
between the S and T each point of S is matched to one or more points of T and vice 
versa [3]. 

Eiter and Mannila [8] reduced the many-to-many matching problem to the minimum- 
weight perfect matching problem in a bipartite graph, and solved it using the Hun- 
garian method in 0(n 3 ) time. The many-to-many matching between two sets on the 
real line is solved by reducing it to the problem of finding the shortest path through a 
directed acyclic graph in 0(n 2 ) time and finally Colannino et al. [T] presented 

an 0(n log n)-time solution for this problem using the dynamic programming method. 

A generalization of the many-to-many matching problem, limited capacity many- 
to-many matching problem, is that in which the number of the points that can be 
matched to each point is limited. Let Cs — {cti, «2, • • • , ot y } and C't = fa, ■ ■ ■ , fa} 
be the capacity sets of the points in S and T, respectively. A minimum-cost limited- 
capacity many-to-many matching is a matching that matches each point s, S S to at 
least 1 and at most on points in T, and each point tj € T to at least 1 and at most j3j 
points in S, for all i,j where 1 < i < y and 1 < j < z, such that sum of the matching 
costs is minimized. Schrijver proved that a minimum-cost limited capacity many- 
to-many matching can be found in strongly polynomial time. A special case of the 
minimum-cost limited capacity many-to-many matching problem is that in which 
both S and T lie on the real line and the cost of matching s, £ S to tj £ T is equal 
to the distance between Si and tj. Panahi and Mohaddes [12] proposed an 0(kn 2 ) 
time algorithm for the one dimensional minimum limited capacity many-to-many 
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Fig. 2.1. The notation and definitions in partitioned point set AU B. 



matching, called ODMLM-matching, where k = max(max\<i< y cti , maxi<j< z (3j ) . 
In this paper we give an 0(n 2 ) time algorithm for this problem. 

The remainder of this paper is organized as follows. Preliminary definitions are in 
section [5] Our algorithm is described in section 01 The matching algorithm presented 
in PP, computing an optimal many-to-many matching between two point sets on the 
real line, is briefly described in section [31 

2. Preliminaries. In this section we describe some preliminary definitions and 
assumptions used in this paper. We represent both the point a and its ^-coordinate in 
the plane using the same symbol a. Let S = {s 1; s 2 , . . . , s y } and T = {t%, t 2l ■ ■ ■ , t z } 
be two sorted sets of points on the real line, that is Si < s% < . . . < s y and ti < t 2 < 
. . . < t z , such that si be the smallest point in S U T. Let S U T be partitioned into 
subsets Aq, Ai, A 2 , . . . such that all points in Ai are smaller than all points in Ai + i 
for all i: the point of highest coordinate in Ai lies to the left of the point of lowest 
coordinate in Ai + \ (Figure l2~Tj) . 

Let A w = {a 1; a 2 , . . . , a s } with ai < a 2 < . . . < a s and A w+ i = {6 l5 b 2 , ■ ■ ■ , b t } 
with bi < b 2 < . . . < b t . We represent \b± — a,| by e,, |6j — 6i| by It is obvious 
that /i = 0. Moreover, ao represents the largest point of A w _\ for w > 0. These 
definitions are presented in Figure HOI In an ODMLM-matching, a point matching 
to its capacity number of points is called a saturated point. 

3. The many-to-many matching algorithm in one dimension . In this 
section, we briefly explain the dynamic programming matching algorithm presented 
by Colannino et al. [1] which determines a minimum-cost many-to-many matching be- 
tween S and T. The time complexity of their algorithm is 0(n log n) for the unsorted 
point sets S and T, but if the point sets are given in sorted order the algorithm runs in 
linear time. Let C{q) denote the cost of a minimum cost many-to-many matching for 
the set of the points {p £ S U T\ p < q}, the algorithm presented in jT| computes C(q) 
for all points q in S U T. In fact, the cost of the minimum many-to-many matching 
between S and T is equal to C(m), where m is the largest point in S U T. In the 
following, we describe the lemmas and the matching algorithm presented in [TJ. The 
proofs of these lemmas can be found in [T] . 

Lemma 3.1. Let a e S,b e T and c G S,d G T such that a < b < c < d. 
A minimum-cost many-to-many matching contains no pairs (a,d), so any matching 
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Fig. 3.1. Suboptimal matchings. (a) (a, d) is a suboptimal matching, (b) (a, b) and (6, d) do not 
both belong to an optimal matching, (c) (a, c) and (b, d) do not both belong to an optimal matching. 

(a, d) in a minimum-cost many-to-many matching, with a < d, satisfies a € Ai and 
d G A i+ x, for some i > fFiaure [3J]a ) J]^. 

Lemma 3.2. Let b < c be two points in S , and a < d be two points in T such that 
a < b < c < d, then both of (a, b) and (b, d) could not be in an optimal many-to-many 
matching, so they are mutually exclusive (Fiaure \3.1b ) fljj. 

Lemma 3.3. Let b < c be two points in S, and a < d be two points in T such 
that a < b < c < d, then a minimum cost many-to-many matching that contains (a, c) 
could not contain (b,d), and vice versa fFiaure \3.1b ) f^. 

Lemma 3.4. In a minimum cost many-to-many matching, each for all i > 
contains a point q i7 such that all points in A^ lying to the left of qi are matched with 
the points in Ai-\ and all points in Ai lying to the right of qi are matched with the 
points in Ai + \ \Tj. 

In fact, we can consider the point qi defined in Lemma 13.41 as a separating point. 
As the separating points are explored, their algorithm easily computes the output. 
They assumed that C(p) = oo for all points p S Aq. Let A w = {ai, 02, . . . ,a s } and 
A w +i = {bi, 62, • • ■ , bt}. Their dynamic programming algorithm's goal is to compute 
C(bi) for each bi 6 A w+ \, assuming that C{p) has been computed for all points p < bi 
in S UT. Depending on the values of w, s, and i five cases arise. 

Case 0: w = 0. In this case two situations arise: i < s and i > s. In the first 
situation the optimal matching is computed by assigning the first s — i elements of 
A to 61 and the remaining i elements pairwise (Figure [3?2k K so C (bi) is computed 
from: 

s i 

In the second situation, i > s, the cost is minimized by matching the first s points 
in Ai pairwise with the points in Aq, and the remaining i — s points in A\ with o s 
(Figure GOb). So : 

s i 

C{bi) = (i- s )e s +Y / e j + J2fr 

Case 1: w > 0, s = t = 1. This case is illustrated in Figure |3~3"r . By Lemma I3T41 b\ 
must be matched with the point a\. So we consider a\ only when it reduces 
the cost of C (61). So: 



C(b 1 )=e 1 +min(C(a ),C(a 1 )). 
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(a) (b) (o) 

Fig. 3.3. (a). Case 1: w > 0,s = t = 1. (b) Case 2: w > 0, s > l,t = 1. (c) Case 3: 
w > 0, s > l,t = 1. 



Case 2: iw > 0, s = 1, £ > 1. By Lemma 13.41 in a minimum cost matching all points 
in A w+ i should be matched with a\ as presented in Figure 13.3b . As case 
1, C (bi) includes C (ai) if a\ covers other points in A w _i; otherwise, C (bi) 
includes C (ao). So the value of C (bi) is: 

i 

C (h) = ]T fi + ^ + min ( c M , C (ax)) . 

Case 3: w > 0, s > 1, t = 1. According to Lemma T3.41 the point ai 6 A w should be 
determined such that all points lying to the left of <Zj are matched to points 
in A w _i and all points lying to the right of ai are matched to points in A w+ \ 
(Figure l3~3b ). So is the point that minimizes the value of C (pi) in the 
equation: 

5 

C(h)= min(Ve J +C(ai_ 1 )) . 

1<2<S ^ ' 

j=i 

Case 4: w > 0, s > 1, t > 1. In this case, we should find the point g that splits to 
the left and right. Let X (bi) be the cost of connecting b\, b%, ■ ■ • , bi to at least 
i + 1 points in A^, (Figure EH^) • Let Sj = J^* =i &j + C(ai-i) for 1 < i < s 
be the cost of connecting the points 04,04+1, . . . ,a s to 6i, plus C(a,_i), and 
Mi = min{5i, S2, ■ ■ • , Si}, then: 

i 

X {hi) = M,-i + fj, !<«<s. 

3=1 



Let y (bi) be the cost of connecting bi, 62, ■ ■ ■ , &j to exactly i points in A, 
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(Figure EOb). then: 

s i 
J=a— i+l j=l 

Finally, let Z (bi) denote the cost of connecting b%, b%, . . . , bi to fewer than i 
points in A w , as depicted in Figure EOb . then we have: 

s i 

Z (bi) = min (V* e h + fj + (i+j — s — X)e s + C (flj-i)) , 1 < L 

s — i-\-2<j<s A — ' z — » 
h=J 3=1 

Note that 

Z (bi) = e s + fi + min(Y (6j_i) , Z 
y (&;) = y (h-i) + e s _ 4+ i + /i + C (a s _ 2 ) - C (a s _ m ) , 

and 

Mj = min(Mi_ x ,Si). 

So the values of X (bi), Y (bi) and Z (bi) can be computed in 0(s + t) time. Note that 

C min(X(fe 4 ),y(^),^(^)) l<i<s 
C(6i)= <^ mm(Y (b s ) , Z (b s )) i = s , 

so the value of C for all 1 < i < t can be computed in 0(s + 1) time. 

4. Our matching algorithm. In this section, we introduce our matching algo- 
rithm which is based on the algorithm of Colannino et al. [T]. It is not hard to see 
that the Lemma 13.31 also applies for an ODMLM-matching. But instead of Lemma 
13.11 Lemma 13.21 and Lemma 13.41 we present other lemmas considering the capacity 
limitations. Let Cap(q) denote the capacity of the point q, i.e. the number of the 
points that can be matched to q. For any point q, C(q,j) denotes the cost of an 
ODMLM-matching for the set of points {p G S U T,with p < q and Cap(q) = j}. 
Note that initially C (pi, k) — oo for all pi 6 S U T and 1 < k < Cap(pi). If m is the 
largest point in S U T, then C(m, Cap(m)) is the cost of an ODMLM-matching. 

Lemma 4.1. Let a G S,b € T and c G S, d G T such that a < b < c < d. If an 
ODMLM-matching contains the pair (a, d), then either b or c is a saturated point. 

Proof. Suppose that the lemma is false. Let M be an ODMLM-matching that 
contains such a pair (a,d), and neither b nor c is a saturated point (Figure l3.1b ). 
Remove the pair (a,d) from M and add the pairs (a, b) and (c,d). The result M' is 
still an ODMLM-matching which has a smaller cost, a contradiction. □ 

Lemma 4.2. Let b < c be two points in S, and a < d be two points in T such 
that a < b < c < d. If an ODMLM-matching contains both of (a,b) and (b,d), then 
the point c is a saturated point. 

Proof. Suppose that the lemma is false. Let M be an ODMLM-matching that 
contains both (a, b) and (b, d), and c is not a saturated point (Figure [3. lb ). Then we 
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(c) 



Fig. 3.4. Case 4-' «i > 0, s > l,t > 1. (a) Computing X(bi). (b) Computing Y(b{). (c) 
Computing Z(bj). 




a b d c 



Fig. 4.1. Suboptimal matchings. (a, c) and (b,d) do not both belong to an optimal matching. 

can remove the pair (b,d) from M and add the pair (c,d): the result M' is still an 
ODMLM- matching which has a smaller cost, a contradiction. □ 

Lemma 4.3. Let a be a point in S , and b, c be two points in T such that a < b < c. 
If an ODMLM-matching contains (a, c), then the point b is saturated with some points 
d in S with d < b. 

Proof. Let M be an ODMLM-matching that contains (a, c), and the point b is 
matched to some points d in S. Suppose that the lemma is false. For d two cases arise: 
b < d < c, and c < d. Assume that b < d < c (Figure |4~T|) . Replace (a, c) and (6, d) 
in M by the two pairs (a, b) and (d,c): the result M' is still an ODMLM-matching 
that has a smaller cost, a contradiction. Now assume c < d, then M contains both of 
(a, c) and (b,d) (Figure 137Tb ) . This contradicts Lemma I3~3l □ 

Lemma 4.4. In an ODMLM-matching, consider Ai sets for alii > 0. Then, each 
Ai contains a point qi, called here the optimal point, such that all points in Ai less 
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a 1 ... °- s — k ... a s 



a l ■■■ a k + l ■■■ as 



"t-s °t-s+l 



Fig. 4.2. Illuustration of Lemma \4--6\ (a) a s < 61. (b) bt < 01. 

than qi are matched to some points b with b < qi, all points in Ai greater than qi are 
matched to some points b' with b' > qi, and qi is matched to either some points b < qi 
or some points b' > q^. 

Proof. If Ai contains a single point, the lemma is true. Now, assume that \A%\ > 1 
and the lemma is false. Let a < b be two points in Ai, the point a be paired with c, 
and the point b be paired with d where b < c and d < a. This contradicts Lemma 
1531 □ 

Lemma 4.5. Let A — {a±, 0,2, ■ ■ ■ , a s } and B = {b±, hi, . . . , bt} be two sets of 
points on the real line, such that the number of the points of B is equal to sum of the 
capacities of the points of A, i.e. t — an< ^ there is an ODMLM-matching 

between A and B. Then \ A\ < \B\. 

Proof. Assume that the lemma is not true and s > t, then since the capacity 
of each point aj e A is ctj > 1, so t < a j- This contradicts the assumption 

Lemma 4.6. Let A — {a\, ai, . . . , a s } and B = {b±, 62, • ■ • , bt} be two sorted sets 
of limited capacity points on the real line ,such that either a s < 61 or b t < a±. Let 
Ca = {<xi, ol%, . . . , a s } and Cb = {Pi, P2, ■ ■ ■ , Pt} be the capacity sets of A and B, 
respectively, such that t < Y]j—i ctj and t > s. Then, we can compute an ODMLM- 
matching between A and B in 0(s + t) time. 

Proof. Either a s < b\ or ai > bt- Obviously, the minimum cost is obtained by 
matching each point in B with as closets point as possible in A. So, the two sets are 
matched as following. 

Assume t < J2j=i a j an d < &i fFigure l4~2k ). Initially, we match the first s 
points in B pairwise with the points in A. This guarantee that each point of A is 
matched at least once. Then, starting from & s +i, we match each point of B to the 
closets point of A, here the largest point of A, that has not been saturated yet. So 
the points in A are matched with the points in B such that the last k points of A are 
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(a) 



(b) 



h l b 2 b 3 b 4 65 



Fig. 4.3. The minimum-cost matching obtained by Lemma \4- 6\ for Ca = {2, 2, 2}: (a) 03 < 61. 
(b) b 5 < "1 • 

saturated, the first s — k — 1 points of A are single matched, and the (s — k) th point 
of A may be neither saturated nor single matched, where k is the greatest integer that 
s — k — 1 + J2j= s -k+i a j < *> i- e - X^= s -/c+i Otf<t — s + fc+1. In fact, a s -k is the 
border between the single matched points in A and the saturated points, such that 
all points in A lying to the left of a s -k are single matched and all points in A lying 
to the right of it are saturated. Figure 14.3a illustrates an example of this situation. 

Now assume t < Xw=i a j and a\ >bp In this situation, that is shown in Figure 
??,we assign the last s elements in B pairwise with the points in A to guarantee 
that each point of A is matched to at least one point of B. Then starting from 
bt- s , we match each point of the remaining t — s points in B, that is 61, ... , 6t_ s , 
to the smallest point in A that has not been saturated. Hence, the first k points 
in A are saturated, the last s — k — 1 points are single matched, and the (k + 1) th 
point may be neither single matched nor saturated where k is the greatest integer 
that s — k — 1 + ^2j = i otj < t, i.e. a j < t — s + k + 1. Figure l4~2b illustrates an 

example of this case of Lemma 14.61 □ 

Theorem 4.7. Let S andT be two sets of points on the real line with \S\+\T\ = n. 
Then, a minimum- cost many-to-many matching with limited capacity between S and 
T can be determined in 0{n 2 ) time. 

Proof. We can compute an ODMLM-matching, by finding the optimal points 
defined in Lemma 14.41 The optimal point of each set A4 for all i separates the points 
lying to the left of the optimal point from the points lying to the right of it. Con- 
sider A w = {ai,a 2 , . . . ,a s } and A w+1 = {61,62, • ■ • , h}- Let C w = {ai,(X2, a s } 
and C w +i = {/3i, 02, ■ ■ ■ , f3t} be the capacity sets of the points in A w and A w+ i, 
respectively . 

Assume that we have computed C (p, k) for all points {p G S U T, with p < 
a s and 1 < k < Cap(p)}, and now we want to compute C(6^,fc) for each 6j G 
A w+ \and 1 < k < Pi. Our algorithm consists of two steps, the primary step and 
the main step explained in the following. 

Primary step, starting from a\, we investigate each aj G A w one by one in the 
ascending order. For each aj G A w that the demands of the points aj,aj+\, . . . ,a s 
can not be satisfied by the points in A w +i, i.e. X)l=i Pk < ( s — j + 1 ) , we seek the 
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Fig. 4.4. ^4 partitioned point set that is used for illustration of the situation where ^^L—i @ k ^ 
(s-j + 1). 



next sets to find enough capacities according to Lemma 14.31 Let C = Ylj—i Pji and 
A = {a,j, aj+i, . . . , aj+c-i} be the C consecutive elements of j4„, starting from ciy 
(Figure |4.4[) . We match etj, dj+i, ■ ■ • , flj+c-i with the points of A w+ \ using Lemma 
14.61 so the last Pi points in A are matched to the point b\, the last P2 points of the 
remaining C — pi points in A are matched to the point b 2 , and so on. Let C be the 
cost of matching the points in A to the points in A w+ i using Lemma |4.6[ then: 

j+c-i t 

C = ^ e fc + ^ Pkfk + min(C (flj-i, aj-i) , C (oj, aj - 1)). 
fe=j fc=i 

Let 

A = {oj+C)Oj+c+i, . . . ,o a } U A w+2 = {a'i,a' 2 , ■ ■ ■ ,a' s >}, 
and A w+3 — {b'i,b' 2 , . . . ,&'{'}• Let C w+ z — {P[,P' 2 , ■ ■ ■ >P'f} be the capacities of the 

points of A w+3 , and C = Y?j=iP'j- Let e 'i = _ a 'il and f'i = \ b 'i ~ If 
A w+ 3 can satisfy the demands of the points in {a'i, a'2, • ■ • , a' s >}, that is C > s', the 
process is stopped. Otherwise, if C < s', we match the first C points in A, that is 
A' = {a'i,a'2, ■ • ■ , a'c}, with the points in A w+3 using Lemma [4.61 So the last 
points of A' are matched to the point b\, the last /3' 2 points of the remaining C — f3\ 
points in A' are matched to the point b 2l and so on. Then, let 

i=i j'=i 

Then, we compare sum of the capacities of the points in A„,+5 with the number of 
the points in {a'c+i, a'c+2, ■ • ■ , a' s '} U ^4to+4, so forth. 

This forward process is followed to find the first set, called Ay , which can satisfy 
the demands of the remaining unmatched points in 

j'-l 

A = {aj,a j+1 , . . . ,a s } U (J A k . 

k=w+2 

Note that it is possible that there does not exist a set such Ay. If Ay exists, we 
should match the unmatched points in A with the points in Ay . 

In fact, by Lemma 14.31 the points in A% can be matched to the points in Ay with 
k < k! for all w < k < j' — 1 (Figure |4~5|) . Figure FO)l shows an example for this 
situation. 

Let 

A = {a"i,a" 2 , ■ ■ • ,a"s"} = {aj,a j+ i, . . . ,a s } U [J A k , 

k=w+2 
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Fig. 4.5. Illustration of the situation Pk < ( s ~ J + !)■ The sets that can be matched 

an ODMLM-matching are connected with a line. 
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Fig. 4.6. An example for the situation y\ /3 fe < (s — j + 1) urat/i Cb = {2, 3, 1, 2, 4}. 



and A,v = {6" l7 fe" 2 ,...,6V}- Consider C A = {a"i, a" 2 , . . . , a'V} and = 
/3" 2 , ■ ■ ■ , P"t"} as the capacity sets of A and A,-<, respectively. Let e"i = 
\b" 1 -a" i \andf" i = \b" i -b" 1 \. 

Depending on the values of s" and t" four cases arise: 

1. s " = l,t" = 1. In this case a'\ is paired with b'\. Let C = C + e"i, then 
we have: 

C(6"i,fc) =min(C(6"i,fc), C") for all 1 < k < f3" v 

2. s" = > 1. In this situation we compute the cost of matching the points 
6"i, 6"2, ■•, b"i' to o"i, and then update the costs of C (b"i',k) for all 1 < i' < 
t": 



C (6'V, k) = mm(C (6'V, k) , C + ^ /"^ + i'e"i) /or a/Zl < fc < /3'V 

3=1 

3. s" > l,t" = 1. In this situation the points a'\, a" 2, ■ ■ ■ , a" S " are matched 
with b'\, let C = C + YfLi then: 



C(b"i,fc) = min(C(&"i,fc), C' + ^Vj) for all s" < k < /3" \. 

i=i 

4. s" > l,t" > 1. In this case we compute the cost of matching a"i,a"2, ■ • ■ , a" s " 
to b"i, b"2, . . . , b"i> for all 1 < i' < t" using Lemma l4~6l We discuss two cases; 
s" < i' or s" > i'. 

First assume s" < i'. In this case either J2j=i a " 3 — or Sj=i a "j > ^ 
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J2j=i a "j = then by Lemma [4J3 we have: 

C' = C' + j2a" j e".+±f" j . 

3=1 3=1 

// 

And if J2i=i a " j > then let A;' be the greatest integer that 



s" - jfc' - 1 + E a"j < 

j=s"-k'+l 



and 



+ E a", \=i'-s" + k'+l- a "i> 

j=s"-k' + l J j=s"—k'+l 

and 

C' = C'+ a" j e" j +me" sll _ kl+ £ <". ■ 

3=a"-fc'+l 3=1 3=1 

according to Lemma l4~6l Then, for all s" <i'< t" : 

C (b"i, , k) = min(C (b" v , k) , C) for all 1 < k < P" v . 

Now assume s" > i'. In this situation either Ylj=i P" j = s " or Y?j=i 0"j > 
s". If Ej=i /3", = then: 

C(&'V, ( 0'V)=™n(C(6'V, ( 8'V), C' + E^ + E^V".) • 

3=1 3=1 

And if X}j=i /^"j > s "; then let fc' be the greatest integer that 

fc' 

3 = 1 

k' \ fc' 

f — fc' — i +E^"i = s " - 1 -E^'i' 

3=1 / 3 = 1 

and 

C" = C" + Ee"3+E/3V"3+^/Vi+ E /'V 

3 = 1 3=1 j = fc'+2 

Then, C (b"i>,k) = mm(C {V' v , k) , C) , 

k = P\, if i' = k' 

for { m<k< /3'V if i' = k' + 1 
1 < k < $" v if k! + 1< i' 
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Main step. Now we start to compute C k) for all 1 < i < t and 1 < k < fa. 
For computing C(bi, k) for each b t S A w+ iand 1 < k < fa two general cases arise: 
Case A: i < a j- I R this case where the demands of the points in A w+ \ can be 

satisfied with the points in A w , we discuss five cases. 

Case A.O: w = 0. This case is similar to Case of [1]. If i = a j> then 

by Lemma 14.61 for all 1 < j < fa we have 

s i 

c {bid) = ^ctkSk + ^fk- 

k=l k=l 

For i < Ylj=i a j we discuss two cases: s < i, or s > i. Assume first that 
s < i. In this situation, we match the points in Aq with the points in 
A\ by Lemma 14.61 So the first s points in A±, that is b\, b%, . . . , b s , are 
matched pairwise with the points in Aq. 

If s < i, starting from 6 s+ i each point of the remaining i — s points in 
Ax is matched to the largest point in A that is not saturated. 
Hence, the last k points of A are saturated, the first s — k — 1 points 
are single matched, and (s — k) th point of A may be neither saturated 
nor single matched, where k is the greatest integer that s — k — 1 + 
Ei= s -fc+i a t <i or J2Us-k+i m <i- s + k + 1. Let 

m = i — I s — k — 1 + ai ] — i — s + k + 1 — on 

\ l=s-k+l ) l=s-k+l 

be the number of the points of A\ that are matched to a s _fc. Then the 
value of C (bi,j) for all 1 < j < fa is: 

s s—k—1 i 

C(bi,j)= ^2 aiei+me s - k + ^ e ' + ft 

i=s-fe+l 1=1 1 = 1 

Now assume that s > i. We identify three cases in this situation: 

1. Y^]=i @j < s - I R this situation, we define C(pi,j) = oo for all 
1 < j < fa- 

2. Y^j=i fij — s - I R this situation by Lemma |4~51 the cost of the match- 
ing is minimized when the last fa points of Aq are matched to the 
point b\, the last fa points of the remaining s — fa points in Aq are 
matched to the point b%, and so on. Let R' = J2j=i e j +S}=i Pjfji 
then: 



oo 1 < j < fa 

R' j - fa 



3. J2]=i Pj > s - By Lemma l4!6l in this situation the optimal matching 
is computed by assigning the last i elements in Aq, a s _i+i, . . . , a s , 
pairwise with the points in A\. Then, starting from a s „i we match 
each point of the remaining unmatched points in Aq to the closest 
point, the smallest point, of Ai that has not been saturated. Let 

s k i 
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k be the greatest integer that 



A- 



i - k - 1 + y^/?j < s, 

3=1 



and 



k \ k 



m 



ft =S-i + k + l-J20: 
3=1 / 3=1 




j = 13., if i = k 

m j — ft if i = k + 1 . 
1 < j < ft if k + Ki 

Case A.l: id > 0, ft = 1, i = 1. In this case, Lemma [4.41 implies that b\ 
must be paired with a s . As Case 1 of [1], we identify two case: (?) a s is 
paired with both b\ and some other points p < b±, and (ii) a s is paired 
with only b±. We choose the matching of minimum cost: 

C (h,k) = min(C (h,k) , e\ + min (C (a s _i, a s _i) , C {a s ,a s - 1)) ) 
for all 1 < k < ft. 

Case A. 2: w > 0, s = 1, i > 1. By Lemma [4.41 all points in A^+i are 
assigned to a s . As Case A.l, other points p < b\ may be assigned to a s . 
So: 

i 

C (bi, k) = min(C (6 i7 k) , fj + iei + min (C (o , a ) ,C (ai,ai - i))) 

3 = 1 

for all 1 < fc < ft. 

Case A. 3: iu > 0, s > 1, ft > 1, i = 1. In this case we should examine all 
points dj G A w and determine the point defined in Lemma |4~41 Let 

s 

C fc = min (S^ei + min(C(ah-i,a h -i),C(a h ,a h -l)), 

s — h+l<h<s * — ' 
z— /i 

then we have: 

C(6i,fc) = min(C(bi,k),C k ) forall 1 < k <ft. 

Case A. 4: to > 0, s > 1, i > 1. In this situation, by Lemma T4. 41 we should 
find the optimal point of A w . 

Let and Z(bi) represent the cost of connecting 61,62, •■•,&» to 

exactly i points, and fewer than i points in A w , respectively. Let X(bi,j) 
represents the cost of connecting 61, 62, . . . , 6, to at least i + 1 points, 
where the number of the points that can be matched to bi is limited to j, 
with 1 < j < (3i. The values of X(bi,j), Y(bi) and Z(bi) are computed 
by Lemma 14.61 as following. 
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Let S[ be the cost of connecting the points a i; a i+1 , . . . , a s to b\ plus 

min (C (aj_i, , C (dj, «j — 1)) for 1 < i < s, so 

s 

Si = y^ej +min(C (aj_i,aj_i) ,(7(0,,^ - 1)) . 
j=» 

Let 

i=l i=fc+2 

be the cost of matching a i; di+i, . . . ,a s to i>i, 62, . . . , bj with j < s — i + 1, 
plus mm (C (a,_i, aj_i) , C (a,, a, — 1)), where fc is the greatest integer 
that 

k 

j-k-l + ^PtfiKs-i + l, 
1=1 

and 



rn 



= s-i+i-j+k+i-j2fofi- 



Let M'ij = min{S'i tj , S'i-ij, S' itj }. Note that if 

i-i 

i=l 

then 

S'ij = 00 for 1 < k < fa. 

Now we compute the values of X (bi,j), Y and Z(6j) for 1 < i < 
min(s, t) and 1 < j < as following: 

X(6 i ,j) = M , a _ i , i , 

and 

s— ii Ot s — i) , C {d s — j+i, Q! s — i+l 1)). 

For the value of Z (bj) we have: 



Z (bi) = min (R ih ) , 

s — i+2<h<s 



where 



R + C (ah-i,ah-i) if k = h 

Rm={ R + min(C(a h -i,ah-i),C(a h ,a h -m} ifk-l = h 
R + min(C(ah-i,a h -i),C(ah,a h -l) if k - 1 > h 
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That in which R is the cost of connecting ah, ■ ■ ■ , a s to b\, . . . , bi, so 

s s—k—1 i 

r= a 3 e i + me s-k + e 3 + E ft> 

j=s— k+l j=h j=l 

where k is the greatest integer for which 

s 

s — h + 1 — k — 1+ aij < i, 

j=s-k+l 



m = i — s + h + k— a j . 

j=s-k+l 

In fact, k = h implies that ah is a saturated point, k — 1 — h implies that 
m points are matched to ah, and k — 1 > h implies that ah is matched 
to a single point. 

Let aj be the smallest point in A w that for which 



Ypi + k>s-j + i, 



then the value of C (bi, k) for all 1 < k < Pi is: 

f min (X (b t , k) , Y (6,) , Z (M , C (b i: k)) ks-j + 1 
C (b h k) = I min (Y (b t ) , Z (h) , C (b t , k)) i = s-j + l 

[ min(minj< h <s (Rih) , C (bi, k)) i > s - j + 1 

Case B: w > 0, i > Ylj=i a j- I R this case, the sum of the capacities of the points 
in A w is less than the number of the points in A w+ \, so by Lemma 14.31 
we should seek the previous sets to find enough capacities for b\ , &2 , • • • , bi 
(Figure |4~7|) . Let C = J2j=i a j> tnen the last C elements in A w+ i, that is 

bi—c+i, ■ ■ ■ i &i i) b%, are matched to the points in A w according to Lemma 

14.61 The cost of this matching is: 



3 = 1 j=i-C+l 

Let 

A w ^ U {h,b 2 , • • • , b^ c } = {b'i,b' 2 , b' v ] with b\<b' 2 < ■■■< b' v , 
A w -2 = {a'i, a' 2 , . . . , a' s '} with a'i < a' 2 < ■ ■ ■ < a' a >, 
C w - 2 = {a'i, a'2, . . • , a' s /}, 

and C = Yf 3 =i a ' 3- 

If sum of the capacities of the points in A w ^ 2 is greater than or equal to 
the number of the points in A w -\ U {bi, b 2 , . . . , bi-c}, that is t' < C , the 
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A w _ 2 








h i ••• b i-c b i-c+i •■■ H 



Fig. 4.7. A partitioned point set used for illustration of Case B: i > ^2°j—i ay. 

process is stopped. Otherwise, the last C elements in {b'i, b'2, ■ ■ ■ , b't>}, that 
is {b't'-c+ii ■ ■ ■ 1 b't'-i, b'f}, are matched to the points in A w ^2 according to 
Lemma WM Let e'i = \b\ — a',| and f\ = \b'i — b'i\, then: 

c = c+ , £^ i + E f'r 

3 = 1 j=t'-C'+l 

Now we compare the sum of the capacities of the points in A w _4 with the 
number of the points in A w s U {b'%, b'2, ■ ■ ■ , &V-C'}) an d so on. 
This backward process is followed until finding the first set, called A w i, that 
can satisfy the demands of the remaining unmatched points in the set 

A W '±i U . . . U A w _x U A w+ i, 

i.e. the sum of the capacities of the points in A w i is greater than or equal to 
the number of the unmatched points in the set A^z+iU. . . U A w -\ U A w+ \. It 
is possible that we do not reach such a set A w i , in this situation C (6j, k) = 00 
for all 1 < k < fii. 

Let A = {b"i, b" 2 , ■ • ■ , &"<'} be the set of the remaining points in Ufc=J,'+i ^k, 
A w i = {a"x,a"2, . . . , a" S "}, a-o the largest point in A w >-±. 
Let C A = {a" 1, a" 2,..., a" s „} and C B = {P"i,f3" 2 , ■■■ ,P"i>} be the capacity 
sets of the points in A w i and A, respectively. Let C (b" y , k') be the cost of 
matching the points in A w > and A using Case A.O, Case A.l, Case A. 2, Case 
A. 3, or Case A. 4, depending on the values of w' , s" , and i' . Then, we have 

C (bu k) = x< mjn, C {b" v , k') + C for all 1 < k < ft. 

Figure illustrates an example of Case B, where A = {si}, A\ = {ti}, 
A 2 = {s 2 ,s 3 },A 3 = {t 2l hM,hMM,C s = {5, 2, 3}, and C T = {2,3,3,1,2,3, 
First, the last 5 points in A 3 are mapped to the points in A 2 , then we seek 
the previous sets to find enough capacities, so we match t 2 with s\. 
Note that 

F(6 i )=y(6 i _ 1 ) + e s _ i+ i + /i+ 
min(C (a s _i,a s _j) ,C (a s - i+ i,a s - i+1 - 1))- 
min(C (a s - i+ i,a s - i+ i) , C (a s _ 4+2 , a s -i+2 ~ 1)), 
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*1 *2 *3 *4 *5 *6 *7 

Fig. 4.8. An example for Case B where t = 7, s = 3 and C$ = {3, 2, 3}. 

and 

-Ri/i = R(i-i)h + .A + efe, 

where is the largest point in A w that is not saturated in the matching corresponding 
to R(i-\)h- And finally, the values of X (bi,j) for all 1 < j < fa is computed in 0(s+t) 
time, since 

S'i-i = S'i + e i + min (C (aj_ 2 , aj_ 2 ) , C (oj-i, - 1)) - 



min (C (aj_i,aj_i) ,C (a,, a, - 1)) , 

and S'i.j = S'ij-i + fj — bk, where bk is the largest point of A w+ \ in the matching 
corresponding to S'ij-i that its degree is more than one. So the values of X (h,j), 
Y (bi), and Z (bi) can be computed in 0(s + t) time. Therefore, we can compute 
C (b i} k) for all 1 < i < t and 1 < k < fa in 0{n) time, and our algorithm computes 
an ODMLM-matching between two point sets with total cardinality n in 0(n 2 ) time. 
□ 

5. Concluding Remarks. The limited capacity many-to-many point matching 
is a many-to-many matching where each point has a limited capacity. In this paper, we 
study the one dimensional limited capacity many-to-many matching problem, called 
ODMLM-matching problem, in which we match two point sets that lie on the real 
line. We provide an algorithm that determines an ODMLM-matching between two 
point sets with total cardinality n in 0(n 2 ) time. The two-dimensional version of this 
problem is open which may be solved using the geometric information. 
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